Casimir Momentum of a Chiral Molecule in a Magnetic Field 
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In a classical description, a neutral, polarizable object acquires a kinetic momentum when exposed 
to crossed electric and magnetic fields. In the presence of only a magnetic field no such momentum 
exists classically, although it is symmetry-allowed for an object lacking mirror symmetry. We per- 
form a full QED calculation to show that the quantum vacuum coupled to a chiral molecule provides 
a kinetic "Casimir" momentum directed along the magnetic field, proportional to its rotatory power 
and the fine structure constant. 



The notion that electromagnetic quantum vacuum 
fluctuations contain energy is widely accepted. When 
coupled to matter they give rise to Van-der-Waal forces 
between neutral objects and Lamb-shifts of energy levels 
in atoms [l|. The attraction between metallic plates is 
undoubtedly the most convincing proof for the existence 
of Casimir energy [^. The empty quantum vacuum is 
Lorentz-invariant with a diverging energy density with 
frequency spectrum w^ and zero momentum. The regu- 
larization of this divergence is still an issue in the context 
of the cosmological constant problem [3|. In all other 
cases the UV catastrophe is not essential when binding 
energies and forces are calculated in a dielectric medium 

a- 

The possibility of the quantum vacuum to contribute 
to the momentum of matter was first raised in Ref. [5| . 
For a macroscopic object exposed to crossed electric and 
magnetic fields, Eq and Bq respectively, a QED correc- 
tion was obtained in addition to the classical Abraham 
momentum density, (e— l)(Eo ABq), which emerges from 
Maxwell's equations and is already controversial by itself 
as concerning the prefactor. Being an observable quan- 
tity, the UV divergence of the momentum obtained in 
Ref. [5| poses a real problem if only because the end result 
depends crucially on an heuristic cut-off. In Refs. [y, |7| 
two of us have shown that the divergence can be handled 
in a full quantum theory that goes beyond the dipole 
approximation. Divergencies were seen to disappear by 
mass regularization and a finite but small Casimir mo- 
mentum remained that is a factor a^ smaller than the 
classical effect. In the quantum theory it is also clear that 
two contributions exist to the (pseudo) momentum. One 
is the familiar sum ^j^ ^ /ik(aj^jake + l/2) associated with 
photons with momentum Kk, and often referred to as 
the "transverse" component Isl]. The second more subtle 
"longitudinal" contribution J2i1i-^i^i) stems from the 
electromagnetic gauge field coupled to electric charges 
{qi} and was found to dominate the Casimir momentum. 
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In this work we apply the same formalism 
Refs. [a, lZ| to address the Casimir momentum of a quan- 
tum object subject to only a magnetic field. The latter 
being a pseudo-vector, symmetry imposes the object to 



be chiral, i.e. we anticipate the relation P^'*'' = ^eBo 
with g some pseudo-scalar quantifying the broken mir- 
ror symmetry of the quantum object. Important is that 
no classical contribution has been reported, so that the 
proposed contribution from the quantum vacuum would 
constitute 100% of the effect. Our theory is nonrelativis- 
tic, and gives a finite magneto-chiral Casimir momentum. 
We find g to be proportional to the product of the static 
rotatory power of the molecule and the fine structure 
constant. 

We propose the simplest molecular model that exhibits 
all necessary features to leading order in perturbation 
theory: broken mirror symmetry, Zeeman splitting of en- 
ergy levels and coupling to the quantum vacuum, and 
we neglect relativistic effects. In our model the optical 
activity of the molecule is determined by a single chro- 
mophoric electron within a chiral object which is further 
simplified to be a two-particle system in which the chro- 
mophoric electron of charge qe = —e and mass mg is 
bound to a nucleus of effective charge qN = e and mass 
rriN 3> me- The binding interaction is modeled by a har- 
monic oscillator potential, V^^ ~ ^{ujlx^+ujyy'^+ujlz'^), 
to which wc add a term Vc = C xyz to break the mirror 
symmetry. The coordinates x, y, z are those of the rel- 
ative position vector, r = r^v — fe, and n = "'-^™'' with 
M ~ 771 jv + TTif,. The center of mass position vector is 
R = [ruMYN + meVe)/M. Vc was firstly introduced by 
Condon et al. [9|,|l0[ to explain rotatory power with a sin- 
gle oscillator model. The parameter C is a pseudoscalar 
which derives from the Coulomb interaction of the sys- 
tem with its chiral environment. When an external uni- 
form and constant magnetic field Bq is applied the total 
Hamiltonian of the system reads, H = Hq + Hem + W, 



withi/o= V -^[p, ~ q,Ao{r,)]^ + V"" + Vc, (1) 
Hem - J2 ^k(aL«k. + h + i^ f d'rBg, (2) 
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In the electromagnetic (EM) vector potential we have 
separated the contribution of the external classical field, 
Ao(ri) = iBoArj, from the one of the quantum field op- 
erator, A(ri). In the sum, al^^ and ate are the creation 
and annihilation operators of photons with momentum 
Kk, frequency cjk = ck and polarization vector e re- 
spectively. The magnetostatic energy of the integral is 
a constant irrelevant to us that we will discard. The sys- 
tem obeying the above Hamiltonian posseses a conserved 
pseudo-momentum [y, |8| , 
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which satisfies [H, K] — 0. Its eigenvalues are therefore 
good quantum numbers. K differs from both the conju- 
gate total momentum, P = Pe -I-P7v, and the kinetic mo- 
mentum, Pkin = MR = P-§BoAr- e[A(rAr) - A(re)] . 
The terms in K can be also arranged as, 



K 



+ PAbr + P 
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where, beside the kinetic momentum, PAbr ~ eBg A r is 
the Abraham momentum and we define the Casimir mo- 
mentum, pCas _ pCas _(_ pCas ^ ^^ ^^ic momentum of the 

vacuum field. P^^s jg composed of a longitudinal Casimir 
momentum, P^^^ — e[A(rjv) — A(re)], and a transverse 
Casimir momentum, P^as _ X^k c ^''('^L'^kc + 5) M- 
Our aim is to evaluate P^^as ^^ ^j^^ ground state of our 
system. The Casimir momentum, we will show, is equiv- 
alent in magnitude to the kinetic momentum acquired by 
the chiral oscillator. 

In the first stage of the calculation we neglect the vac- 
uum field. The conserved pseudo-momentum of the re- 
sultant Hamiltonian, Hq, is Q = Pkin + eBg A r, with 
[iJo,Q] = and continuous eigenvalues Q. The uni- 
tary operator U= exp [i(Q — |Bo A r) • R/H] maps the 
Hamiltonian Hq into Hq =\J^Ho\J, which conveniently 
separates internal and external motion, 

^" " ^^' + ^^^ + ^""^ + ^c + Vz + 0{Bl). (6) 

In this equation p is the conjugate momentum of r, p = 

fi{pN/mN - Pe/rrie), and Vz = 2^(r A p) • Bo is the 
Zeeman potential with u* = "'"™'' . 

The ground state of the Hamiltonian Hq is, up to order 
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y cyclic permutations. 



Correspondingly, the ground state of Hq is |f7o) =U|f^o), 
with a pseudo-momentum Qo to be fixed. The fact that 
(rio|r|f^o) — implies also that Qo is the kinetic momen- 
tum of the oscillator in absence of vacuum. In the above 
equation the states \nxnynz) refer to the eigenstates of 
the harmonic oscillator Hamiltonian. The dimensionless 



parameters are, Bq 



eS' 



4ti* ^/uTjZJj: 



with i =/= j J^ k and i ^ k, 



C = g^^ „u = "'-"j Here the 

indices i,j, k take on the three spatial directions, x, y, z. 
The 77 factors are assumed to be small quantities which 
quantify the anisotropy of the oscillator. They all have 
to be nonzero for the optical activity of the molecule to 
survive rotational averaging. In the following, all our 
calculations restrict to the lowest order in C, Bq and if^. 

Next, we couple the oscillator to the vacuum field and 
evaluate the Casimir momentum, (P^as^ _ (JJ7|pCas|J-J^) 
Here, \^) is the ground state of the oscillator when cou- 
pled to the vacuum field. It will be computed applying 
up to second order perturbation theory to |ilo) with the 
interaction potential W . In general. (P^as^ ^jjj depend 
on the pseudo-momentum of |rio), Qq. In the presence of 
the vacuum, Qo is not the kinetic momentum of the sys- 
tem but the quantum number of the operator K instead. 
It can be shown that (P) = Qo - fBo A (r) - (Pj''"). 
By inserting this relation into Eq.(j4]) we find, (K) = Qg. 
For simplicity we take Qo = 0, which corresponds to 
an oscillator that would be at rest in absence of vac- 
uum, and cvluate the expectation value of the opera- 
tors in Eq.([5]). The induced dipole moment, (er), van- 
ishes, and so docs (PAbr)- The rest of the terms yield, 
(Pkin) + (Plf^*^) 4- (P5^^) = 0, for Qo = 0. This 
means that a nonzero value of the Casimir momentum 
implies a nonzero kinetic momentum in the molecule en- 
tirely due to its interaction with the vacuum field. The 
transverse momentum, (P^^"^) = ^j, ^ hk.{n\a]^^a\^f\Vl) ^ is 
found to be finite but small. The longitudinal momen- 
tum, (P?''*') = e(Sl|A(rAr)-A(re)|ri), gives the dominant 
contribution that we will evaluate in detail first. 

At ©(CBo) and lowest order in the coupling constant 
e, it suffices to compute |fi) applying first order pertur- 
bation theory to \^o) . We use the U-transformcd states 
and the U-transformed potential. 
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to arrive at, 
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where AA = A(rjv) — A(re) and Eq ~ h{ujx+^y+^z)/2. 
The intermediate atomic states, /, are eigcnstatcs of Hq 
and may have a priori any pseudo- momentum Q. Zeros 
in the denominator arc avoided in the summation. It 
turns out that only the terms proportional to p • A in VF 
yield a nonvanishing contribution to (Pp'**'). Writing the 



EM quantum field in Eq.® as usual [12 
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with V a generic volume, c the speed of light and eg the 
vacuum permittivity, passing the sums over Q and k to 
continuum integrals and by summing over polarization 
states we arrive at. 
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where I is the unit dyadic and the kinetic energy 
fi,^fc^/2M in the denominator is a consequence of the 
momentum recoil of photon to atom, Q = — k. When 
the oscillatory exponentials on the left of the Hamilto- 
nian operator are moved to the right, two types of term 
appear. To the first type belong those terms in which 
the products of the oscillatory exponentials cancel. They 
correspond to processes in which a photon is created and 
annihilated at the same particle (i.e., electron or nucleus). 
The second class contains terms in which the product of 
the exponentials amount to e**''''. They correspond to 
processes in which a photon is created and annihilated 
at different particles. The first kind of terms arc the 
dominant ones. The second kind of terms are a factor 
yj EqI [ic? smaller and will be neglected. 
The dominant terms are. 
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where —[rrie O m^v] means that the same expression 
within the square brackets must be evaluated exchanging 
me with mjq and subtracted. H^'^ =p'^/2fi + V^'-' is the 
harmonic ocillator Hamiltonian. Note that, since both 
the chiral interaction Vc and the Zeeman effect Vz feau- 
turing in Hq are treated as perturbations to the harmonic 
oscillator, the denominator in Eg. ljlOp is expanded up to 
order VcVz in Eq. ([TT|) . Hence, only terms at 0{CBo) 
must be retained. When moving the oscillatory expo- 
nentials from the left to the right of the Hamiltonian 
operator in Eq. (|10p , in application of canonical commu- 
tation relations the momentum p in Hq gets shifted, 
p — >■ p ± "^^j'° hk. As a result, the kinetic energy in 
the denominators of Ea. ([TT|) becomes h^k'^ /2me^N and 
an additional term ±/ik • p/yTig ^y shows up there. The 
latter is a Doppler shift term which has been neglected 
inEq.([Tll)- 

In Eg. pTj) . the factor which accompanies (Jlnjplpo) is 
a UV divergent integral. It was shown in Refs. [a, l7[ that 
it can be interpreted as a mass renormalization factor. 
In our case this divergence does not have any effect on 
the Casimir momentum since (f7o|p|0o) = 0. All the 
other integrals in Eq. (|lip contain more than one factor 
quadratic in k in their denominators, which makes them 
fully convergent. When integrating in A: up to cx3 the 
leading contribution of them all is logarithmic. The end 
result can be expressed as, 
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where Eijk is the three-dimensional Levi-Civita tensor. 
Note that r/*''-' = —if''. We thus need the anisotropy 



of the oscillator quantified by rj^'^ to make the Casimir 
momentum nonzero. The momentum is still not averaged 
over the random orientation of the chiral object. To do 
so we must project (P^as^j qj^^q ^]^g magnetic field vector 
and extract the trace, 
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where the subscript rot denotes rotational average. The 
logarithmic dependence on Eas. (jl2ll3p may be disputed 
since it results from wave numbers larger than miy^ec/h. 
It is possible that relativistic corrections are needed, 
which lie outside the scope of the present Letter. 

The transverse Casimir momentum, (P^^'^), corre- 
sponds to the expectation value of the operators in the 
sum over modes of Eq.([5]). Its evaluation involves a sec- 
ond order pertubation calculation whose details will be 
published somewhere else ll| . The end result is approx- 
imately, 

,pca.^ _ -0.98Ce3Bo 

\-t^_L /ro 



XAA'n'^ceQm'lijJxi^yUJ 
;p|P^^)™t/ln(mAr/me). 
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For instance, if we take for niN the mass of a car- 
bon atom, {P'j^^)rot is an order of magnitude less than 
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Eq.(fT3|) and Eq.(fT4|) give simple formulas for (P 

and {P^^)rot in terms of the chiral parameter C, the 
magnetic field and the different oscillator frequencies. We 
will now express them in terms of the static optical ro- 
tatory power and electric polarizability of the molecule, 
which can be easily calculated for this simple model and 
which are both directly observable. To this aim we derive 
equations for the electric and magnetic dipole moments 
induced in the state |f2o) of the molecule by external 
electric and magnetic fields of frequency w, E^ and B^ 
respectively. For this we apply first order perturbation 
theory to |r2o) using the potential W in which the quan- 
tum field must be substituted by a classical EM vector 
potential from which the external electric and magnetic 
fields derive. After performing the rotational average we 
can write for the molecule's electric and magnetic dipole 
moments, respectively, 

{d)rot = aE^u+l'Bo^^B^-|3^B^, (15) 

{rxv)rot = aA/B^ + CBoAB^+7BoAE<^-t-/3E^, 

where we have adopted Condon's notation Q in SI units. 
In this equation ue and um are the electric and mag- 
netic polarizabilities respectively. C, describes the mag- 
netic Faraday effect, jS is the molecular rotatory factor 
responsible for the natural optical activity and 7 gives 
rise to the magneto-chiral anisotropy. In agreement with 



Ref. [13j they all behave as ■ 



for w ;» Ei^/h whereas 



they tend to constant values in the static limit, a; — > 0, 
where we find. 
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In this formula I? is a rational polynomial of the natural 
frequencies of the oscillator which approximates 2? « 1 
for not too large anisotropy. By comparing the above 
expression with Eas. (|13ll4p we conclude that the total 
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Casimir momentum, (P ^^)rot = (P| irot 
can be written approximately as. 
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[ln(mAr/me) + l]eBo, (17) 



where a is the fine structure constant. Thus, we 
confirm our conjecture, p^as _ geBo, with g = 
~i^ a (0) P'-^ {^^N I'Tde) -f 1]. We speculate that, apart 
from constants of order unity, this expression is model- 
independent. 

The presence of h in Eq. (fT6)) deserves a comment. The 
fact that /3 is proportional to fi is a consequence of the 
quantum formalism. For a given set of natural frequen- 
cies, all the atomic lengths in the problem are determined 
by quantum mechanics. In particular, /3(0)/a£;(0) is a 



length, a fraction of the electronic Compton wavelength, 
that we identify with a chiral length, l^h- Therefore, we 
can write P '^^^ ^ a^lch'^Cj where wp = eBo//i is the 
cyclotron frequency. 

Next, we compare our result with earlier work on 
Casimir momentum wich is based on the application of 
the fluctuation-dissipation theorem and worked out in 
the dipole approximation 14|, |l5|. For our model that 
approach would yield. 
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This expression is analogous to the formula obtained 
in Refs. [5|, ll6| for the transverse Casimir momentum 
transferred to a homogeneous magnetoelectric medium. 
In contrast to the quantum calculation, the integral of 
Ea. ([T5)) diverges quadratically for k^ax -^ oo and needs 
a UV cut-off to yield a finite result. Upon comparing the 
result of this integral with the quantum result of Eq. (fT4)) 
we have identified this cut-off with the inverse Compton 
wavelength, kmax ~ ^- The details of the calculation 
will be published elsewhere |11| . 

Finally, we estimate the orders of magnitude of the 
Casimir momentum transferred to a single molecule of 
an actual chiral compound. To this aim we derive values 
for /3(0) and as(0) from experimental parameters. Let 
us consider a molecular medium with a number density 
p. At leading order in p and far from resonances the ro- 
tatory power of the medium, (/?, relates to the molecular 
rotatory factor /3 as, (/? = 47r^p/3(0)/A^eo, and the refrac- 
tive index is ri = 1 -I- pa£;(0)/2eo. As an example we take 
the chiral compound 2-octanol, CgHigO. According to 
Table I of Ref. ^], at A = 4800A, n == 1.43 and its spe- 
cific rotatory power is [(/j]''^*°°=15.46 deg/dm/(g/cm^). 
Its mass density is 0.82 g/cm'^ and its molecular mass 
is 130.2 u. We find, l3{0) = 3 x lO-'^^T&dC^m^ /3, 
aEiO) = 2 X lO-^^m^C^/J. Considering that the chiral 
center is one of the carbon atoms and using Bq = lOT, 
we find (P^'''')rot = 1.4x 10"^^kgm/s. This corresponds 
to an average velocity of 0.6 nm/s approximately, which 
lies within the scope of current experimental observations 

In this Letter we have carried out a quantum com- 
putation of the Casimir momentum transfcrcd from the 
quantum vacuum to a chiral oscillator subject to an ex- 
ternal magnetic field. The result is free of divergences, 
linear in Bq and proportional to a, which proves the gen- 
uine QED nature of the effect. Our calculation is based 
on a simplified molecular model which allows us to ex- 
press the final result in terms of observable parameters, 
i.e., the molecular rotatory factor and the electric polar- 
izability. Ea. (|T7)) is the main result of this Letter. A 
quantitative prediction has been made whose verification 
is accessible to current experiments. 
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